A Method to Modify RMT using Short-Time Behavior in Chaotic Systems 
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We discuss a modification to Random Matrix Theory eigenstate statistics, that systematically 
takes into account the non-universal short-time behavior of chaotic systems. The method avoids 
diagonalization of the Hamiltonian, instead requiring only a knowledge of short-time dynamics for a 
chaotic system or ensemble of similar systems. Standard Random Matrix Theory and semiclassical 
predictions are recovered in the limits of zero Ehrenfest time and infinite Heisenberg time, respec- 
tively. As examples, we discuss wave function autocorrelations and cross-correlations, and show 
how the approach leads to a significant improvement in accuracy for simple chaotic systems where 
comparison can be made with brute-force diagonalization. 

PACS numbers: 05.45.Mt, 03.65.Sq 



The statistical structure of chaotic wave functions has 
been a key topic of investigation from the early history 
of quantum chaos and wave chaos physics, and its study 
is essential for improved understanding of resonances, 
transport, and long-time dynamics in non-integrable sys- 
tems [l|. Random Matrix Theory (RMT) 0, which of- 
ten serves as an adequate zeroth-order approximation for 
wave function statistics in the absence of integrability, 
describes a statistical ensemble of Hamiltonians having 
no preferred basis. Within RMT, eigenstates are sim- 
ply random vectors either in the full Hilbert space or 
in the subspace given by energy and other conservation 
laws. For a quantum particle in a slowly-varying poten- 
tial, a wave function then behaves locally like a random 
superposition of plane waves of fixed wave number, as 
discussed by Berry Q . 

As a universal theory, RMT specifically excludes 
system-specific behavior associated with dynamics, 
boundary conditions, or interactions. Well-recognized 
deviations from random wave function statistics are asso- 
ciated with boundary effects [J, @ , finite system size [H , 
unstable periodic orbits @, diffusion and two-body 
random interactions in many-body systems 0, Q • Much 
progress has been made in understanding such deviations 
in various situations of physical interest, for example 
chaotic wave function correlations in Husimi space associ- 
ated with classical dynamics [9j and realistic mesoscopic 
S-matrices arising from a simple diffusive ray picture of 
wave propagation (lpj. In particular, semiclassical meth- 
ods [111 ] have proven very successful in quantifying the 
effects on wave functions of boundaries [j, [f| and peri- 
odic orbit scars Q. However, the limit implied by semi- 
classical approximations may not always be achievable 
or relevant in describing actual experiments. For exam- 
ple, an analysis of electron interaction matrix elements 
in ballistic quantum dots shows that even for thousands 
of electrons in the dot, several statistical quantities of 
interest typically exceed random wave predictions by a 
factor of 3 or more; for other quantities the random wave 
model fails even to predict the correct sign [ijj (see also 
0)- 

In some situations, e.g., [12j . brute force diagonaliza- 
tion of the Hamiltonian may be used to obtain correct 



statistics for the stationary or long-time behavior, but 
for very large Hilbert spaces, such as those that arise 
in many-body situations, diagonalization is likely to be 
impractical. Even where it "works" , diagonalization is 
unlikely to produce much intuition about the relevant 
physics, and must be repeated for each new Hamilto- 
nian. In fact, individual eigenstates of a chaotic Hamil- 
tonian are highly sensitive to perturbations of the system, 
particularly for multi-particle systems. The statistics of 
such systems are far more robust and remain accurate 
for small perturbations. 

Our goal here is to present a system and basis- 
independent way of supplementing RMT with short time 
dynamical information, that eliminates the need for diag- 
onalization of the Hamiltonian, and that provides greatly 
improved accuracy over RMT and semiclassical methods 
for finite systems with a finite Ehrenfest time. 

To enable direct comparison with RMT, let us consider- 
fully chaotic (ballistic or diffusive) dynamics without 
symmetry on an iV-dimensional Hilbert space with eigen- 
states To avoid ambiguities in the definition of we 
assume a non-degenerate spectrum. Typical quantities of 
interest, then, are functions of the amplitudes (a|£) for 
any physically- motivated basis state \a), which may be 
a position or momentum state, a Slater determinant, or 
more generally an eigenstate of some zeroth-order Hamil- 
tonian. With the normalization X^^Li l( a l£)| 2 = 1j the 
simplest and first non-trivial moment of these amplitudes 
is given by the local inverse participation number (IPR) , 
which measures the degree of localization at |o): 
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(1) 

varying from P aa = 1 in the case of perfect ergodicity 
to P aa = N for perfect localization. For two arbitrary 
states we have 
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P ab = Nj2\(a\0\ 2 \(b\0\ 2 = Nlim w dt\(a\b(t))\ 



(2) 
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Obviously, higher-order moments and in general the en- 
tire joint distribution of the eigenstate intensities may be 
considered (e.g., [Hj]). We may also relax the requirement 
that only pure states such as |a)(a| act as probes, and in- 
stead measure the structure of chaotic eigenstates using 
any desired self-adjoint operator a [lH ]. Operator probes 
(of phase space size greater than or smaller than h) will, 
for example, be particularly helpful in the study of hi- 
erarchical eigenstates in a mixed chaotic-regular phase 
space 16]. Again, without loss of generality we may 
adopt the normalization Tr a = 1 . Eq. @ becomes 
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(3) 

with the autocorrelation P aa as an obvious special case. 

In the semiclassical limit N — > oo, averages of the 
form (jTJ)- ([3]) may be obtained using short-time dynam- 
ics; specifically for discrete-time dynamics we have 



where 



pab r^^pab^ 5 



P ab (t) = \(a\b(t))\' + (aHt)){b(t)\b), 



(4) 



(5) 



where \a) and \b) may be any two states (identical, over- 
lapping, or orthogonal) Here and in the following, 
7TT indicates an ensemble average. If desired, the ensem- 
ble may be selected so that all realizations possess the 
same short-time dynamics P ab (t), in which case the aver- 
age on the right hand side of ((U) is superfluous. The cut- 
off time t must be long compared to the ballistic or diffu- 
sive Thouless time (so as to include all the non-universal 
dynamics), and short compared to the Heisenberg time, 
which scales with N. No distinction is made in f4} be- 
tween non-universal short-time revivals that indicate de- 
viations from RMT in the eigenstate statistics and the 
0(1/ N) short-time revivals that are present already in 
the context of RMT. As a result, ^ systematically over- 
estimates corrections to RMT, and violates probability 
conservation ^2 h P ab = 1 given a complete basis \b) for 
any r > 0, with the violations growing linearly as t/N. 

We now notice that the problematic aspects of (01 for 
finite system size ./V can be eliminated by introducing a 
t- and (a|6)-dependent prefactor: 



P ab « C< o|b> (r) / dtP ab {t) 



(6) 



where in particular C^'^(t) = N/At converges to the 

exact answer as t — + oo. To fix cffl b \ we apply RMT to 
Eq. ^ and obtain 



pab ^ pab 
r ~ RMT " 



(7) 



Eq. ([7]), and its natural extensions to higher-order mo- 
ments (e.g., (P af >)™) and operator expectation values 
(e.g., P a P) are a key result of this paper. Stationary 
eigenstate properties of an quantum chaotic system or 
ensemble of systems may be fully described by a com- 
bination of short-time dynamics for that system or en- 
semble, in combination with exact results from RMT, 
without any need for matrix diagonalization. Reassur- 
ingly, Eq. ([7]) yields exact results in three limits of in- 
terest: (i) the RMT limit where P ab (t) = P^ MT {t) and 
thus P ab = -PrmT' ( u ) the semiclassical limit N/t — > oo, 
where we recover (|4]), and (iii) the limit where an infinite 
amount of dynamical data is available as input, t — > oo. 
More importantly, as we will see in the examples below, 
Eq. ([7]) and its extensions provide reliable approxima- 
tions to exact diagonalization in situations far from any 
such limit, i.e., for finite-size systems far from universal- 
ity, and where the only input is short-time dynamics on 
the scale of a Lyapunov time. 

Short-time overlaps P ab (t) needed as input to Eq. ([7]) 
may sometimes be known analytically, as in the case of 
periodic orbit scars, while in more general situations the 
short-time dynamics for a given system of interest is eas- 
ily obtainable numerically, to any desired time scale r. 
The RMT factors in ([7]) and its generalizations may be 
treated entirely analytically. For example, for arbitrary 
| a) and \b) we have standard results in the absence of 
time reversal symmetry (GUE or CUE) 
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(8) 



while for a general self-adjoint operator a we obtain 



paa 
RMT 



\ 1 / 



where A4 are the eigenvalues of a Ai — 1). 

Similarly, RMT dynamical overlaps may be expressed 
exactly using RMT eigenstate statistics and the RMT 
spectral form factor, e.g., 



P&T® = J^PrIt + £(e^'-^<)RMT 



x (\(a\0\ 2 \(bW)\ 2 + (a|c>(£|o}(aa)(&|e>) (10) 

V / RMT 

For discrete-time dynamics, described by the CUE en- 
semble, which will be relevant for the numerical examples 
below, we have 



1 



for t = 



2 

N+l 



PrmtW = (1 + l<#>| 2 ) x { M for 1 <\t\ < N , 

for \t\ > N 

(11) 

and analogous results for self-adjoint operators are ob- 
tained by spectral decomposition, as in @. 
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Thus, eigenstate statistics for a chaotic system or en- 
semble of systems may be unambiguously obtained with- 
out diagonalization, as in (|7|), by combining exact RMT 
results with easily obtainable short-time dynamical in- 
formation for the system or ensemble of interest. 



of the statistics, we have added the functions t(p) and 
v(q), which are random within a small region near the 
edges of the phase space (\p\ > 1/2 — A and \q\ > 1/2 — A) 
and zero elsewhere. In the following, we set A = 0.1, but 
the results have no significant dependence on A. 
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FIG. 1: The inverse participation ratio P aa for a Gaussian 
distribution centered on a short periodic orbit with instability 
exponent Ao = 0.5 is computed by direct diagonalization and 
compared with the short-time dynamical prediction given by 
Eq. 0. Here the system size is N = 32 and the Gaussian 
distribution has size s = 0.5 (Left panel) or s = 0.25 (Right 
panel) . Convergence to the exact result is observed when the 
dynamical calculation includes information about times r up 
to 2 or 3 in units of the local Lyapunov exponent Ao. The 
RMT value Prmt = (1 + s~ 1 )iV/(iV + 1) and the semiclassi- 
cal result Pgg = (1 + s- 1 )(7V/(iV + 1)) j^Z-oo sech(A i) are 
shown for comparison. All quantities appearing here and in 
subsequent figures are dimensionless. 

We now discuss a few illustrative examples, using as 
our model the paradigmatic example of a quantized pe- 
riodically kicked Hamiltonian [TtJ 



H(q,p,t)=T(p) + V(q) £ *(* 



(12) 



on the compact phase space (q,p) G [— 1/2, 1/2) 2 . The 
kinetic and potential terms are chosen to produce a fully 
chaotic map (perturbed cat map [Lsj ) 



Tip) 



K 
K 



cos(2-7rp) + t(p) 



V(q) = -~q 2 - C0s(2^<j) + v{q) 



(13) 
(14) 



where the parameters to and K control the chaoticity 
of the system: the dynamics is fully chaotic for m > 
| if | and the instability exponent of the shortest periodic 
orbit at q = p = is A = cosh -1 (l + (to - K) 2 /2) w 
to — K for m — K <C 1. To break time reversal and 
parity symmetries, and also allow for ensemble averaging 
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FIG. 2: The inverse participation ratio for a pure Gaussian 
wavepacket (s = 1) is computed exactly and compared with 
the dynamical prediction of Eq. J7| using dynamical informa- 
tion up to times r = 2AQ -1 and 4AQ -1 , where Ao = 0.25 is the 
local Lyapunov exponent. Results are shown for various val- 
ues of the system size N. The semiclassical and RMT limits 
are also shown for comparison (see Fig. [T] caption). 

We begin by considering the inverse participation ra- 
tio P QQ ; where a is the Weyl transform of a Gaussian 

— 2 / 2 _ 2/2 

distribution p(q,p) ~ e q l(J i p > a p centered on the pe- 
riodic orbit. We define s = <r q a p /h. Then in the special 
case s = 1, a is a projection onto a minimum uncertainty 
Gaussian wave packet, while more generally a represents 
a mixed initial state. Typical results are shown in Fig. [TJ 
where the dynamical prediction of Eq. (0 for several val- 
ues of the cutoff time r is compared with exact values 
obtained by brute-force diagonalization. We note that 
the dynamical prediction begins at the RMT limit for 
t = 0, as it must, and quickly converges to the exact 
stationary answer at 2 or 3 Lyapunov times. Fig. [5] il- 
lustrates the relationship between the exact value of the 
inverse participation ratio, the dynamical prediction, and 
the limiting RMT and semiclassical approximations, as 
the system size N is varied. Here we note significant 
deviations from the semiclassical answer even when N 
takes values of 100 or greater; these deviations are well 
reproduced in the dynamical calculation. 

As another example, we consider wave function inten- 
sity correlations P ab for position states \a), \b). Since 

-4j ^2ab=i P ab = 1 is given by wave function normal- 
ization when t — > Qo , we focus on the first interesting 
moment, the variance 



w =^T,( pab y 2 - 1 



(15) 



a . b 



W is a simple measure of non-uniformity in infinite- 
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time transport [19(, and ranges from W = for per- 
fect crgodicity to W = N — 1 for prefect localiza- 
tion. We note also that interchanging the roles of eigen- 
states and basis states, W may be equivalently writ- 
ten as the variance of the interaction matrix elements 
P«' = NJ2a=i \( a \0\ 2 \{aW)\ 2 between eigenstates |£) 

and i.e., W = j^Ef.e' (-P"') 2 - L The statistics 
of such interaction matrix elements in chaotic systems 
frequently appear in applications ranging from quantum 
dot conductance in the Coulomb blockade regime [12| to 
controlling directional emission properties in microcavity 
lasers Hi- 
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FIG. 3: The interaction matrix element variance W is com- 
puted exactly (Eq. (|T5j)) and compared with the short time 
prediction with tA = 1. Here A = 0.125. The RMT 

result and the semiclassical limit Wsc — 1/N are also shown 
for comparison. 



We note here that the intensity correlators P ab pre- 
dicted by Eq. ([7]) are not guaranteed to satisfy the nor- 
malization condition p^f tl i ,ail = 1, that holds for 
the exact correlators. This normalization is only guar- 
anteed for (Eq. ([7])) when r — > oo. In order to pre- 
dict W we already need knowledge of the short-time dy- 
namics, P ab (t), for every pair of initial and final states 
| a), 1 6); therefore with little added computational ef- 
fort we may achieve exact normalization and further 
improve the convergence with r, simply by rescaling 



A - 



Fig. [3] shows that the semiclassical and RMT predic- 
tions are very similar for the system we consider here, 
and both deviate significantly from the exact results for 
finite N. Our method, including rescaling, converges to- 
ward the exact answer very quickly, on the order of the 
Lyapunov time, even where the RMT prediction is off by 
a factor of 2 or 3. The accuracy can be improved further 
by calculating the dynamics for longer times r. 

We have developed a method that improves on RMT 
eigenstate statistics for chaotic systems by systemati- 
cally incorporating short-time dynamics. The method 
is conceptually appealing, computationally simpler than 
brute-force diagonalization, and significantly more accu- 
rate than RMT or the semiclassical limit for realistic 
systems. The approach can be easily extended to con- 
sider symmetry effects (including time reversal symme- 
try), mixed phase space [2l| . and resonance wave function 
statistics in open systems. 



We again combine short time dynamics and RMT to 
calculate the variance of the interaction matrix elements, 
similarly to Eq. ([7]), 
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